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Abstract 

We study the abihties of the Fermi-LAT instrument on board of the Fermi mission to simul- 
taneously constrain the Milky Way dark matter density profile and some dark matter particle 
properties, as annihilation cross section, mass and branching ratio into dominant annihilation 
channels. A single dark matter density profile is commonly assumed to determine the capa- 
bilities of gamma-ray experiments to extract dark matter properties or to set limits on them. 
However, our knowledge of the Milky Way halo is far from perfect, and thus in general, the 
obtained results are too optimistic. Here, we study the effect these astrophysical uncertainties 
would have on the determination of dark matter particle properties and conversely, we show 
how gamma-ray searches could also be used to learn about the structure of the Milky Way 
halo, as a complementary tool to other type of observational data that study the gravitational 
effect caused by the presence of dark matter. In addition, we also show how these results would 
improve if external information on the annihilation cross section and on the local dark matter 
density were included and compare our results with the predictions from numerical simulations. 
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1 Introduction 



One of the mysteries of modern physics is the presence of an unknown non-luminous and non- 
baryonic matter in the Universe (see Refs. [IHS] for reviews), the so-cahed dark matter (DM). 
However, although there are many pieces of evidence in favor of the existence of DM and which 
indicate that it constitutes about 80% of the total mass of the Universe [GH^, it is worth noting 
that they each infer DM's presence uniquely through its gravitational influence. In other words, we 
currently have no conclusive evidence for DM's non-gravitational interactions. Thus, the origin and 
most of the properties of DM (mass, spin, couplings) remain unknown. 

Although many particles have been proposed as DM candidates (see, e.g., Refs. |4i|9j for a 
comprehensive list), a weakly interacting massive particle (WIMP), with mass lying from the GeV 
to the TeV scale, is one of the most popular ones. WIMPs can arise in extensions of the Standard 
Model such as supersymmetry (e.g., Ref. [1]), little Higgs (e.g., Ref. [lOj) or extra-dimensions 
models (e.g., Ref. ) and are usually stable and thermally produced in the early Universe with 
an annihilation cross section (times relative velocity) of {av) 3 x 10~^^ cm^/s, which is the 
standard value that provides the observed DM relic density. Thus, WIMPs are cold DM (CDM) 
candidates. 

Within the framework of CDM, as favored by observations, structure forms hierarchically bottom- 
up, with DM collapsing first into small halos, which then accrete matter, merge and eventually give 
rise to larger halos |12fll4) . Within this paradigm, the properties of luminous galaxies are expected 
to be related to those of DM halos [15| . which carry with them gas that cools and collapses to their 
center to form galaxies. 

Over the past decades, the progress in high-resolution N-body simulations has permitted a better 
understanding of the structure of CDM halos. It has been shown that density profiles of CDM halos 
can be reasonably well described by an universal form, independent of the halo mass, which predicts 
cuspy mass distributions. Although the first analytical analysis of this kind predicted that the radial 
distribution of DM follows a simple power-law |16| . it was later shown with N-body simulations that 
the slope of the DM distribution is not constant and changes with the distance to the center of the 
halo. It was Navarro, Frenk and White (NFW) who first proposed a simple two-parameter fitting 
formula to describe spherically averaged DM density profiles |17|ll8j. However, further simulations 
with improved numerical resolution started to show some systematic differences |19H32) . mainly in 
the innermost regions of CDM halos. Thus, in order to reproduce several of these features, as the 
gradual shallowing of the density profile towards the center of the halo, an improved three-parameter 
fitting formula was proposed [33], the so-called Einasto profile |34) . 

Nevertheless, these prescriptions only represent the mean of all simulated halos for a given mass 
at a given redshift. The scatter with respect to these mean values |35H40) is thought to be related 
to the different halo formation histories and to the evolution of the expanding Universe |18 |l41fH9] . 
Thus, these average values do not necessarily describe accurately the DM halo of our own Galaxy. 
Indeed, it is not clear that the Milky Way is a prototypical spiral galaxy. For instance, while the 
presence of the two Magellanic Clouds points to a more massive halo, this would also imply higher 
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chances for a recent merger that could have destroyed the thin galactic disc |50) . 

However, probing the inner structure of the Milky Way DM halo is a very challenging task 
and current data do not allow to distinguish among different density profiles |51|I52). This is 
particularly important for the astrophysical detection of DM, as these searches are critically sensitive 
to the structure of the Milky Way DM halo. Direct detection experiments, searching for signals of 
nuclear recoil of DM scattering off nuclei, depend on the local DM density and velocity distribution. 
On the other hand, indirect searches look for the products of DM annihilation (or decay), which 
include antimatter, neutrinos and photons. In the case of DM annihilation, the luminosity depends 
quadratically on the DM density and the largest uncertainties come from the region where we 
expect the highest (neutrino or gamma-ray) signal, i.e., the galactic center (GC). Hence, a better 
understanding of the Milky Way halo is required to determine the potential reach of current and 
future experiments. Conversely, astrophysical searches, and in particular indirect DM gamma-ray 
searches, could also be used to learn about the distribution of DM in the innermost regions of 
the Milky Way halo, beyond the limit of convergence of numerical simulations. Thus, in case of a 
positive signal, these searches could be useful not only to learn about the DM nature, but also as a 
complementary tool to constrain the Milky Way DM density profile, which may not be necessarily 
representative of the general characteristics obtained in numerical simulations. 

In this work, we study the abilities of the Large Area Telescope (Fermi-hAT) instrument on 
board of the Fermi Gamma-ray Space Telescope (Fermi) mission [53] to constrain the DM density 
profile and the effect the uncertainties in the DM distribution could have to constrain DM properties, 
as annihilation cross section, mass and branching ratio into dominant annihilation channels. We 
consider the potential gamma-ray signal in a squared region with a side of 20° around the GC 
{\l\, \b\ < 10°) and, in order to model the relevant gamma-ray foregrounds, we use the numerical 
code GALPROP [54] and the latest Fermi-L AT observations [551(56]. 

During the last years, different approaches have been proposed to determine the WIMP DM 
properties by using indirect or direct measurements or their combination [571489) . In addition, 
the information that could be obtained from collider experiments would also be of fundamental 
importance to learn about the nature of DM [90fll05) and could also be further constrained when 
combined with direct and indirect detection data [58p 98 |[T06lI108) . Therefore, in order to show how 
well the DM density profile could be determined by the combination of different type of experiments, 
we will also make use of potential information that might be obtained on some parameters by other 
complementary means in addition to the gamma-ray signal. 

The paper is organized as follows. In Section [2] we describe the most relevant (for this work) 
component of the gamma-ray emission from DM annihilation in the GC. In Section [3] we state 
some of the results obtained in numerical N-body simulations regarding the density profiles of CDM 
halos and we also briefly review the current knowledge of some of the most important parameters 
that determine the structure of the Milky Way halo. We describe the ingredients of our analysis in 
Section U] and present the results in Section [5j Finally, we draw our conclusions in Section [6l 
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2 Gamma-rays from DM annihilations 



The differential flux of prompt gamma-rays generated from DM annihilations in the smooth Milky 
Way DM halo after the hadronization, fragmentation and decay of the final states and coming 
from a direction within a solid angle Afi can be written as |112) 
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where Rq is the distance from the Sun to the GC, Pq = p{Rq) is the local DM density, {av) is 
the thermal average of the total annihilation cross section times the relative velocity, rriy- is the 
DM mass, the discrete sum is over all DM annihilation channels, BRj is the branching ratio of 
DM annihilation into the i-th final state and dN^/dE^ is the differential gamma-ray yield of SM 
particles into photons for the i-th channel, that we simulate with the event generator PYTHIA 
6.4 |113| . which automatically includes the so-called final state radiation (photons radiated off the 
external legs). The dimensionless quantity J, which depends crucially on the DM distribution, is 
defined as 

^(^) = 7^17^ / / P{ris^^)fds, (2) 

AiZ UqPq JaQ J\os 

where the spatial integration of the square of the DM density profile p{r) is performed along the 
line of sight within the solid angle of observation AQ. More precisely, the distance from the GC is 
r = ^ Rq — 2s Rq cos tp + s'^, and the upper limit of integration is Smax = i^tir ~ ^ ^q) ~^ 
Rq cosV', where ^ is the angle between the direction of the GC and that of observation. The virial 
radius of the Milky Way DM halo is defined as R^h = [3 Mvir/(47r A^ir Pcrit)]^^^) where Mvir is the 
virial mass of the DM halo, pcrit is the critical density of the Universe and Avir = ISvr^ + 82 (fim — 
1) - 39 (Qm - 1)^ is the virial overdensity [HI] with fim = 0.258 [7J the non-relativistic matter 
contribution to the total density of the Universe. It is also customary to define the concentration 
parameter as Cvir = Rvh/'>'s, where is the scale radius. 

3 Milky Way DM density profile 

During the past two decades, the improvement in numerical simulations has led to the discovery of 
a universal internal structure for spherically averaged DM halos |17yi8| that can be well described 
by a density profile (NFW parametrization) with two (free) parameters given by 

[r/rs) [I + {r/rs)\^ 

However, further simulations with improved numerical resolution started to show some systematic 
differences, mainly in the central regions of CDM halos, that have created a vivid debate [191432) . 
The recent Via Lactea II simulations |110j seem to partly verify earlier results, with a NFW density 



^Throughout this work we neglect the contribution due to substructure in the Milky Way halo, which could 
enhance the gamma-ray flux from DM annihilation by a factor of <10 [109H111"| . 
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profile. On the other hand, the Aquarius project simulations seem to favor a slightly different 
parametrization |33 |I115M118) . which is less cuspy towards the center of the Galaxy than the NFW 
profile and has three (free) parameters, the so-called Einasto profile |34) . 



which is the same relation between slope and radius that defines Sersic's model |119| . but applied 
to space density instead of the projected surface density of galaxies [120J. 

In order to observationally determine the parameters of the DM density distribution, different 
constraints are considered. For instance, several dynamical tracers have been used to determine the 



modeling of the radial density and velocity dispersion profiles of distant halo tracer populations, as 
satellite galaxies, globular clusters and horizontal branch halo stars on far orbits. However, while 
the mass of external galaxies can be determined with reasonable precision, the mass of the Milky 
Way remains uncertain within a factor of ~ 6 (5 x 10^^ M© < Mmw < 3 x lO^^ j^g) [mtilM ISTl 
[S2 tll29l4l32| due to the lack of observational data for luminous tracers in its outer regions where 
DM dominates. Within the innermost ~ 50 — 60 kpc of the Milky Way, on the high mass end we 
find Mmw(< 50 kpc) = Mmw(50) = (5.4l°;^) x 10^^ Mq [122], whereas a lower mass was recently 
obtained by estimating the circular velocity curve up to 60 kpc, Mdh(< 60 kpc) = Mdh(60) = 
(4.0 ± 0.7) x 10^^ M0 [T28|. Here Muwir) and Mdh(r) refer to the total mass (dark and visible) of 
the Milky Way and the DM halo mass within a radius r, respectively. Below we will compare these 
estimates with the masses for the profiles we consider, but we stress that the former estimate refers 
to the total mass within a given radius and not only to the mass of the DM halo. The galactic 
disc and bulge (visible matter) are estimated to contribute with a total mass about an order of 
magnitude lower, ~ (5 — 7) x 10^^ Mq |133pi34|ll32j . which would bring the former estimate closer 
to the latter one. 

The local DM density (pq) is another important parameter in the description of the DM halo. 
This is, as well, a crucial parameter in direct DM detection experiments and in the search for 
neutrinos produced by DM annihilations at the center of the Sun or the Earth. Commonly, the way 
to estimate the local DM density follows the approach of Ref. |135) and it is based in constructing 
three-component (disc, bulge and dark halo) models for the galaxy and confront them against 
observational data. It is usually assumed to be = 0.3 GeV/cm^, although approximately with a 
factor of two of uncertainty. However, this estimate is just a rough central value within the range 
obtained by several analyses from old dynamical constraints and using cored DM profiles [136) . 
Moreover, it is curious to note that in these early estimates, to the best of our knowledge, there is 
no justification for using that particular value, which is not even an unweighted average, but on the 
contrary, Pq was usually quoted to be slightly larger, e.g., Pq = O.Sllg^^ GeV/cm^ [136] . 

Recently, a few new analyses have been performed using new data. Some of them follow the 
approach of Ref. [135) but with a large set of new observational constraints of the Milky Way [511 

■^More precisely, these analyses estimate the mass up to some radius, typically tens of kpc, and then infer the total 
mass from a model of the Milky Way. 




(4) 



Milky Way mass □, and most estimates have 



come from escape velocity arguments or from Jeans 
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[52 t l79 |[T371ll32| and in Ref. [138] a new technique is proposed such that only local observables are 
used and there is no need of a model for the galaxy. In general, they tend to obtain slightly larger 
values than the canonical one, />© = 0.39 ± 0.03 GeV/cm^ [M], = 0.32 ± 0.07 GeV/cm^ [79], 
P0 = 0.43 ± 0.15 GeV/cm^ [138], p© = 0.40 ± 0.04 GeV/cm=' [l32]. Let us also note that DM halos 
are not predicted to be perfectly spherical and this possible flattening could imply a larger DM local 
density |136pi39] . In this work we will take /j© = 0.4 GeV/cm^ as our default value. 

Correlated with the use of the DM local density is the knowledge of the distance to the GC (Rq), 
which is also an important parameter to determine the structure of the Milky Way. It is usually 
assumed that Rq = 8.5 kpc. However, this is based on an old recommendation by the International 
Astronomical Union, after the value for the unweighted average of Rq obtained using different 
estimates between 1974 and 1986 |140] . Nevertheless, after a more careful statistical analysis of 
these results, which tried to account for statistical and systematic errors and the covariance among 
the different methods, a lower value of Rq = 8.0 it 0.5 kpc was proposed almost twenty years 
ago |141) . Many new measurements were performed during the following decade, with an average 
value Rq = 7.80 it 0.33 kpc |142j . consistent within errors with previous estimates. 

During the last years, and using new astrometric and radial velocity measurements of stellar 
orbits around the massive black hole in the GC, new analyses seem to point to a slightly larger 
value, Rq = 8.4 ± 0.4 kpc [US], Rq = 8.33 ± 0.35 kpc [HI], Rq = 8.2 ± 0.5 kpc [SS]. In addi- 
tion, based on trigonometric parallax observations of star forming regions, a similar result has been 
derived, Rq = 8.4 it 0.6 kpc |146) (see however Ref. |147| ). All in all, even with these seemingly 
converging results, recent measurements still span in the range between Rq = 7.2 it 0.3 kpc |148) 
and Rq = 8.7^q^ kpc [149]. In this work, we will take Rq = 8.3 kpc as our default value. 

From Eqs. ([3|) and ([3]), we see that the DM density profile is determined by two parameters in the 
NFW case {pq, Vg) and three for the Einasto profile (pq, Vg, a). Nevertheless, one could equivalently 
use any set of (free) parameters to describe the DM density profile with a given parametrization. For 
instance, in N-body simulations, it is usually convenient to use the virial mass and the concentration 
parameter instead of the local DM density and the scale radius. On the other hand, it turns out 
that different analyses, which make use of the results from N-body simulations, have shown that 
the structural properties of DM halos depend on the halo mass |18 1l35I[371l39fl49| . In general, the 
correlations among different parameters also embed the dependence on the cosmological parameters. 
In particular, these analyses have found an anti-correlation between Mvir and Cvir, which is especially 
sensitive to the value of erg (a measurement of the amplitude of the linear power spectrum on the 
scale of 8 Mpc). The median relation between Mvir and Cvir (assuming a NFW profile) at 
redshift z = and for relaxed halos using the cosmological parameters from the 5- year release of 

■^These are halos in dynamical equilibrium, in general with smooth density profiles and reasonably well fitted 
by NFW (and Einasto) profiles. However, the degree of relaxedness could be very important to determine halo 
properties |150) . 
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the WMAP mission [7], is given by [39] 

where h = 0.72 0. Hence, this type of correlation would reduce further the number of parameters to 
describe the DM density profile to one (to two in the analogous case of the Einasto profile). However, 
it is important to bare in mind that Eq. ([5]) represents the mean concentration for a given virial 
mass, so it does not account for the intrinsic scatter about these mean relation [35fl37| l39 l l40) . The 
distribution of the logarithm is approximately Gaussian and at la is typically ux^c^.^^ — 0.3 |35H40| 
(see however Ref. [T51j where instead of a log-normal distribution of halo concentrations, a Gaussian 
distribution is found and where also the uncertainty in DM annihilation from the Milky Way halo 
that arises from the distribution of halo concentrations is discussed). In the analysis that follows, we 
will also consider the possibility that the Milky Way satisfies this mean relation (Eq. ([5])) obtained 
in numerical N-body simulations. 



4 Analysis 

The Fermi telescope [53J was launched in June 2008 for a mission of 5 to 10 years. Fermi-LAT is the 
primary instrument on board of the Fermi mission and it performs an all-sky survey with a field of 
view FoV = 2.4 sr, covering an energy range from below 20 MeV to more than 300 GeV for gamma- 
rays, with an energy and angle- dependent effective area which is approximately A^s ~ 8000 cm^. 
We consider a constant effective area, but we notice that different effects could reduce the exposure 
after a more detailed description of the detector is considered, which however is beyond the scope 
of this paper. Fermi-LAT's equivalent Gaussian la energy resolution is ~ 10% at energies above 
1 GeV. On the other hand, for normal incidence, 68% of the photons at 10 GeV are contained within 
^ 0.25°. In the following analysis, we consider a 5-year mission run, an energy interval from 1 GeV 
to 300 GeV and a square region of observation with a side of 20° around the GC (|/|, \h\ < 10°). 
In our simplified modeling of the detector, we do not explicitly compute the effects of the energy 
resolution and of the point spread function. However, in order to avoid that this affects significantly 
our results, we divide the energy interval into 20 evenly spaced logarithmic bins and consider 10 
concentric 1° bins. Our aim is to illustrate the issues here discussed and serve as a starting analysis 
for more detailed ones. 

In order to study the capabilities of the Fermi-hAT experiment to simultaneously constrain 
some DM particle properties and the Milky Way DM density profile, we begin by modeling the 
background. We consider the three components contributing to the high-energy gamma-ray back- 
ground: the diffuse galactic emission (DGE), the contribution from resolved point sources (PS) 
and the isotropic gamma-ray background (IGRB). We model the DGE with the GALPROP code 
(v54) |54j by means of the WebRun service [152] and we take the conventional model, which re- 
produces reasonably well the recent measurements by Fermi-LAT, at least at intermediate galactic 
latitudes 10° < |6| < 20° and up to 10 GeV [153|I154] . The other two sources of background are 
taken from the Fermi-hAT published results [55 1 I56 | . We refer the reader to Ref. [61] for further 
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details on the gamma-ray foregrounds we consider. However, let us note that Ref. |61| makes use 
of the Fermi Science Tools, both to model the background and the signal. We have checked that 
our simplified modeling of the detector and backgrounds does not introduce significant changes in 
our results and allow us to substantially reduce the time of computation without modifying our 
conclusions. 

During the last years, different approaches have been considered to constrain WIMP DM prop- 
erties by using indirect searches, direct detection measurements, collider information or their combi- 
nation |57fll08| . These measurements are complementary and constitute an important step toward 
identifying the particle nature of DM. However, astrophysical searches are hindered by large un- 
certainties in the structure of the Milky Way DM halo and, although different profiles have been 
considered, it is customary to assume a given DM density profile in order to study the capabilities 
of gamma-ray experiments to extract some DM properties |57fl62) (for attempts to also fit the DM 
density profile see Ref. |57) , but assuming the annihilation channel to be known and Refs [155|I156] , 
but for a very large annihilation cross section). 

In this work, we discuss the Fermi-LAT's abilities to constrain the DM density profile and some 
DM properties after 5 years of data taking. We evaluate the effect the uncertainties in the DM 
density profile could have on the determination of the DM mass, annihilation cross section and the 
annihilation channels. Concerning the DM density profile, when we consider an Einasto profile, we 
fit its three parameters: the DM local density Pq, the scale radius and the index a; whereas for 
the case of a NFW profile we fit Cvir and Mvir- In principle, as with respect to DM properties, the 
analysis should include all possible annihilation channels, but this would render the present analysis 
very much time consuming. Nevertheless, in practice, they are commonly classified as hadronic (soft 
channels) and leptonic channels (hard channels). Thus, for simplicity, when simulating a signal, we 
only consider two possible (generic) channels, t~^t~ or By using this simplification, we can 

reduce the number of total free parameters to six for the Einasto profile and five for the NFW profile: 
the mass, rriy., the annihilation cross section, {(Jv), the branching ratio into channel 1, BRx(2) (or 
equivalently into channel 2, BR2(i) = 1 — BRx(2)) and the parameters describing the DM density 
profile (pqjVs, a) or (Mvir, Cvir) for the Einasto or NFW profile, respectively. In order to focus on the 
determination of these parameters and to avoid adding more parameters to the analysis, we assume 
perfect knowledge of the background. Hence, in this regard, the results here presented should be 
taken as the most optimistic ones. We use the function defined as 

j=l i=l v"statjjj -r {UsysJij 

where Sij{9) represents the simulated signal events in the i-th energy and j-th angular bin for 
each set of the parameters 9 = [pq, Vg, a, m^, {(Jv), BR^(ft)) for the case of the Einasto profile or 



''Note that for these channels, and for the energy range of interest here, the contribution from the electrons 
and positrons produced in DM annihilations to the gamma-ray spectrum via inverse Compton scattering off the 
ambient photon background, is subdominant with respect to prompt gamma-rays [61| . Also note that for the con- 
sidered masses, electroweak radiative corrections to DM annihilations with gauge boson bremsstrahlung do not 
contribute [75711166] . 
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Table 1: Summary of the parameters and relevant information for each of the figures in 
Section{^ We take as default values: Rq = 8.3 kpc, DM annihilation into a pure bb final state and 
{av) = 3 X 10~^^ cm^/s. The Einasto (NFW) profile is used when a value for a is (not) indicated. 
When priors are added, criog^av) 

= 0.15 is always assumed and for Pq two cases are considered: 
(^pQ = 0.04 GeV/cm^ and Up^ = 0.2 GeV/cnr' . All the figures assume 5 years of data taking. 



9 = (Mvir,Cvir, {(^v) ^ BR^(f,)) for the NFW profile, and Sj^{d^) is the assumed observed signal 
events in that bin with parameters 6^ = (^p%, r^, a^, m^, {av)^ , BR^^^^^^ for the Einasto profile or 

9^ = ^M^jj., c^jj., m^, (av)^, BR^^^-j^ for the NFW profile. The statistical errors are given by 

{<7st.t)l = Sj^ + B,, , (7) 

where Bij is the number of background events in the i-th energy and j-th angular bin. On the 
other hand, in the analysis of Ref. [56], the rms of the residual count fraction between Fermi-hAT 
data and the DGE model for energies above 200 MeV is larger than the expected statistical errors 
by a factor of 2.5. Hence, baring in mind that the number of events in each bin is in general 
dominated by the background events, we conservatively take the systematic errors to be three times 
the statistical ones, i.e., agys = 3<Tstat) in order to account for the uncertainties in the determination 
of the background. We expect however that this will improve in the near future. 

5 Results 

In this section, we first consider the Einasto profile and show the results for two different DM masses, 
= 25 GeV and m^. = 80 GeV. We show the capabilities of the Fermi-LAT telescope to constrain 
the three parameters that describe the DM density distribution. In addition, we also show the effect 
our ignorance about the halo profile has on the determination of some DM properties. This could 
be compared with the results obtained in Ref. [61], where a known DM density profile was assumed. 
Next, we also study a NFW profile for = 25 GeV for different values of the parameters of the 
DM profile and plot the results in the (Mvir, Cvir) plane. In this way, one can readily compare our 
results to the relations obtained in numerical N-body simulations, e.g., to Eq. ([5|). Moreover, we 
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[GeV] a 

Figure 1: Fermi L AT abilities to simultaneously constrain DM properties (left panel) 
and the Milky Way DM density distribution (right panel) after 5 years. We consider DM 
annihilation into a pure bb final state (BR^(^iy-j = 0), (crv) = 3 x 10^^^ cm^/s and = 80 GeV. 
For the DM density profile we consider an Einasto profile with pq = 0.4 GeV/cm^, Vg = 20 kpc, 
a = 0.17 and Rq = 8.3 kpc. The black crosses indicate the values of these parameters. In both 
panels, blue (orange) regions represent the 68% CL (90% CL) contours for 2 dof. 

also study the impact of adding some priors on some of the relevant parameters, which could be 
obtained in the future using different data, as {av) and Pq. In Table[T]we summarize the parameters 
used and the relevant information in each of the figures we describe below. 

5.1 Constraining DM properties and the DM density profile 

For the first four figures, Figs.[lHll we consider an Einasto profile (Eq. Q) with our default values for 
the local DM density and the Sun distance to the GC {pq = 0.4 GeV/cm^ and Rq = 8.3 kpc) and we 
further set the other two parameters to = 20 kpc and a = 0.17. These values for the scale radius 
and the index a are the ones commonly assumed in studies of DM indirect detection. Moreover, 
this value of a approximately represents the mean value obtained from N-body simulations for dark 
halos with virial masses of the order of 10^^ Mq |118) . Concerning the DM particle properties, we 
assume DM annihilation into a pure bb final state and an annihilation cross section (times relative 
velocity) (av) = 3 x 10~^^ cm^/s. The left panels of the four figures show the results for pairs of 
DM properties after marginalizing with respect to the rest of the parameters and the right panels 
show the results for the parameters that describe the DM density profile. 

In Fig. [H we depict the Fermi~hAT reconstruction prospects after 5 years. The blue regions 
and the orange regions correspond to the 68% and 90% confidence level (CL) contours for 2 degrees 
of freedom (dof), respectively. However, note that there are very small differences between these 
contours, which is due to a very steep x^- 

From the left panel one can extract two important pieces of information. On one side, unlike 
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Figure 2: Fermi L AT abilities to simultaneously constrain DM properties (left panel) and 
the Milky Way DM density distribution (right panel) after 5 years. Same as Fig. but for 
= 25 GeV. 

what happens if the DM density profile is assumed to be known |61| . in the present case the value of 
the annihilation cross section would remain completely unknown. This can be easily understood by 
the fact that by also fitting the DM density profile, the quantity /Jq J is allowed to vary and along 
with (<Tf), that quantity sets the normalization of the potential signal (see Eq. ([T|)). Thus, there is 
a clear degeneracy between them, which cannot be broken by the energy or angular binning of the 
signal. On the other hand, the capabilities for the reconstruction of the DM mass and annihilation 
channel do not substantially worsen with respect to the case of a known DM density profile (but 
with a single angular bin) |61| . This is clearly due to binning the signal both spatially and in energy, 
as there are more "data" points and the energy spectrum obviously depends on the DM mass and on 
the annihilation channel. Hence, the DM mass could be determined with a ~ 25% uncertainty and 
the branching ratio could simultaneously be constrained to be BR.^({,) < 0.25 at 90% CL (2 dof). 

From the right panel of Fig. [T| one realizes that it could be possible to establish a strong cor- 
relation between and a at a high CL (note again that the \^ is very steep). It is important to 
stress that this could be achieved by only using gamma-ray observations with no further constraints 
on the relevant parameters and basically thanks to the angular binning. There are also correlations 
among these parameters and p©, although they are much weaker. Nevertheless, due to the degen- 
eracy between /Jq J and (erf), it turns out that the actual values of the parameters that describe 
the DM density profile cannot be determined with only gamma-rays. Although being a very weak 
constraint, a could be found to lie in the range a = (0.12 — 0.30) at 90% CL (2 dof). 

In Fig. [2] we show the corresponding results for the case of m-^ = 25 GeV. The same trends as 
in Fig. [1] are observed, although due to larger statistics, slightly better results would be obtained. 
We see that the DM mass could be determined with a ~ 10% uncertainty and simultaneously, the 
branching ratio could be constrained to be BR^(;,) < 0.10 at 90% CL (2 dof). In addition, the index 
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Figure 3: Fermi L AT abilities to simultaneously constrain DM properties (left panel) and 
the Milky Way DM density distribution (right panel) after 5 years, with the addition of 
external information on pQ and (crv) (see text). Same parameters as in Fig.Ul In both panels, 
the dark (light) regions represent the case of Ciog^^^^ = 0.15 and = 0.1/90 ((TpQ = 0.5 Pq). As 
in Fig.m blue (orange) regions represent the 68% CL (90% CL) contours for 2 dof. 



a could be bounded to lie in the range a ~ (0.13 — 0.27) at 90% CL (2 dof). We also note that the 
correlations between parameters get tighter, but the degeneracy between Pq J and (av) cannot be 
broken, which would prevent us from being able to determine neither the annihilation cross section 
nor the exact parameters that describe the DM density profile. 

So far, we have shown results that could be obtained by just using gamma-ray observations. 
However, we already have some information about some of the parameters entering the fit, as the 
local DM density (pq) |511l52p79pi371ll32pi38) . and one would also hope that other parameters might 
be constrained in the future, as the annihilation cross section {{crv)) and the DM mass (m^) |70H108]. 
In Figs. [3] and m we show how our knowledge about the structure of the Milky Way dark halo could 
be largely improved if we add to the gamma-ray analysis, constraints coming from other sources on 
Pq and {av). 

For the local DM density we will consider two different priors. On one side, we follow the 
common and conservative approach of assuming that it is determined within a factor of two and 
set a 50% error, cTpg = 0.5 Pq. On the other hand, we also consider a more optimistic constraint as 
recently obtained in different works |51 |I132) . cXpQ = 0.1 Pq. 

Concerning some of the DM properties, we will only set a prior on the annihilation cross section. 
Although there are good prospects of constraining the DM mass by using direct searches and collider 
experiments, we will not implement this potential constraint in what follows. As shown above, a 
gamma-ray observation could allow a quite good determination of the DM mass if the DM particle 
is relatively light (m^ < 100 GeV). Indeed, we have checked that a further external constraint on 
this parameter would not add substantial information. However, as also discussed above, there is 
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Figure 4: Fermi L AT abilities to simultaneously constrain DM properties (left panel) and 
the Milky Way DM density distribution (right panel) after 5 years, with the addition of 
external information on pQ and (av) (see text). Same as Fig.\^hut for = 25 GeV. 



a degeneracy between /Jq J and (crw), so in order to break it as much as possible, not only is the 
information on important, but a constraint on {av) is also needed. Hence, we assume that the 
annihilation cross section might be determined within an order of magnitude at 3(T with future 
collider data |98] and set a very optimistic prior o"iog(o-i)) = 0.15, where {av) is expressed in cm^/s. 
We note that a weaker constraint on this parameter would not improve significantly the results as 
compared to the case with no priors, so in this regard, information on both /Oq and (erf) is crucial 
to reduce the allowed parameter space to determine the structure of the Milky Way DM halo. 

In order to implement these priors, we add a penalty factor to the for each parameter (with 
Gaussian error), which takes into account the external information. The new reads 



Xprior 



H,i2, (s« (9) - SI," («»))■ 



=1 i=l 



(cTstat)-j + (0-| 



+ 



Pi 



'PQ 



+ 



log{av) — \og{av)^ 



(8) 



where again, in log{av), {av) is expressed in cm^/s. 

In Figs. [3] and U we assume the same parameters as in Figs. [1] and [21 respectively, but now 
we add the external information on /Jq and {av) we have just discussed. In both figures, the blue 
(orange) regions correspond to the 68% CL (90% CL) contours for 2 dof. In these figures, the case 
of our most optimistic set of priors (ciog^o-u) = 0-15 and = 0.1 /9q) is represented by the dark 
regions, whereas the slightly less optimistic case (uiog^o-^,^ = 0.15 and ap^ = 0.5 Pq) is depicted as 
the light regions. 

In both figures, we see that there is no further improvement on the determination of the DM mass 
or the annihilation branching ratio, as gamma-rays do not provide more stringent constraints on {av) 
than what we assume to get from collider experiments or DM direct detection searches. Likewise, the 
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"data" on gamma-rays do not add further information on the local DM density. Nevertheless, for the 
case of niy. = 80 GeV (Fig. [3]), the index a could be constrained to lie in the range a ~ (0.13 — 0.22) 
and a ~ (0.12 — 0.25) at 90% CL (2 dof) and for our most and less optimistic priors, respectively. 
Although modest, this could be an useful piece of information. On the other hand, could be 
only weakly constrained for the case with our most optimistic priors, — (8 — 60) kpc at 90% CL 
(2 dof). 

Similar improvements are obtained for a lighter DM, as shown in Fig. [J] for = 25 GeV. Very 
small differences are found between the results for these two DM masses. The allowed range for 
the index a would turn out to be a ~ (0.14 - 0.21) and a ~ (0.13 - 0.23) at 90% CL (2 dof) and 
for our most and less optimistic priors, respectively. On the other hand, with our most optimistic 
priors, the allowed region for the scale radius would be r<j ~ (8 — 56) kpc at 90% CL (2 dof). 

5.2 Comparison with numerical simulations 

The parametrization commonly used when studying the structural properties of DM halos by means 
of numerical N-body simulations is the NFW profile (Eq. (|3])). Thus, in order to compare our findings 
with their results, we consider this parametrization in the next three figures. 

In Figs. [5H3 we show the results for different values of the parameters that describe the DM halo 
(we again set Rq = 8.3 kpc), that is the local DM density and the scale radius or equivalently, the 
virial mass and the concentration parameter. Indeed, for the sake of comparison with simulations, 
we show the results in the plane (Mvir,Cvir) @, after marginalizing with respect to the three DM 
particle parameters we study. However, we do not show here the prospects to constrain these 
parameters, as the results are very similar to those obtained for the Einasto profile and depicted in 
the left panels of Figs. [IHll In the next figures, we assume the same DM mass, rriy. = 25 GeV and, 
as in Figs. [IHll we consider DM annihilation into a pure bb final state and an annihilation cross 
section (times relative velocity) {av) = 3 x 10~^® cm^/s. 

The left panels of Figs. [SHZlshow the capabilities of the Fermi-LAT experiment to reconstruct 
the halo properties after 5 years. For these panels no external priors are assumed, analogously to 
Figs. [Hand [21 As in the Einasto case, the blue and orange regions correspond to the 68% and 90% 
CL contours for 2 dof, respectively. Again, we see that there are very small differences between 
these contours due to a very steep x^- 

The right panels of Figs. [5HZ1 depict the results in the case of adding external information on Pq 
and {(Jv), as in Figs. [3] and [H Likewise, the blue (orange) regions correspond to the 68% CL (90% 
CL) contours for 2 dof. In these panels, the case of our most optimistic set of priors {cri^g(^av) — 0-15 
and (JpQ = 0.1 Pq) is also represented by the dark regions, whereas we show the less optimistic case 
(a"iog(o-t,) = 0.15 and = 0.5/9©) as the light regions. 

In all three figures, in both panels, the red solid curves represent the mean value and la scatter 
band of the concentration parameter obtained in numerical simulations (see Eq. ([5]) and discussion 
below). The central black dashed line indicates our default value for the local DM density, Pq = 

*We only show the range of Mvir currently favored. 
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Figure 5: Fermi L AT abilities to constrain the Milky Way DM density distribution after 
5 years (keeping correlations with DM particle properties) . We consider DM annihilation 
into a pure bb final state (BR^(ij-^ = Oj, {av) = 3 x 10^^^ cm^/s and = 25 GeV. We assume a NFW 
profile, Pq = 0.3 GeV/cn^ , Rq = 8.3 kpc and the relation in Eq. ([S]) fcvir = 9-2, Vg = 38.6 kpc and 
Mvir = 2.5 X 10^^ j. The black crosses indicate the corresponding values in the (M^ij-, Cvir) plane. 
In both panels, blue (orange) regions represent the 68% CL (90% CL) contours for 2 dof. The left 
panels show the case with no external priors and in the right panels two external priors are adopted: 
'^log(o-u) = 0-15 and = 0.1 /?© for the dark regions and (Tiog^o-^^ = 0.15 and = 0.5 Pq for the 
light contours. Red solid lines indicate the mean concentration parameter obtained from numerical 
simulations \49^ and the la scatter hand ((T\ac^;^ = 0.3j; black dashed lines indicate bottom-up 
Pq = 0.2 GeV/cm^, 0.4 GeV/cm^ and 0.6 GeV/cm^. 

0.4 GeV/cm^, whereas the lower black dashed line indicate Pq = 0.2 GeV/cm^ and the upper one 
Pq = 0.6 GeV/cm3. 

We first consider in Fig. [S] the case of Pq = 0.3 GeV/cm'^ and impose the relation Mvir — c^ir for 
the mean obtained in numerical simulations |49] (see Eq. (I5|) ) . In this way, the two-parameter profile 
gets reduced to a one-parameter parametrization. Thus, a local DM density Pq = 0.3 GeV/cm^ 
and Eq. 1^ imply Mvir = 2.5 x 10^^ and Cvir = 9.2. The scale radius is predicted to be 
= 38.6 kpc, which is much larger than the value of 20 kpc commonly used in the literature of DM 
indirect detection. As we see from Fig. [5l the gamma-ray data alone would not be able to set any 
significant constraint on the DM density profile. Moreover, although adding external information 
on Pq and (av) would substantially reduce the allowed parameter space, yet a large fraction would 
still be unconstrained even in the most optimistic of our cases. In particular, no bounds on the 
virial mass could be set with this technique. 

In general, the larger the scale radius (the smaller the concentration parameter) the lower the 
potential signal from DM annihilations. In Fig. |6]we consider Pq = 0.4 GeV/cm^ and = 20 kpc, 

^Let us note that if pq — 0.4 GeV/cm'^, then Eq. ([Sjl imphes Mvir = 6.1 x 10^^ M©, which is probably too large. 
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Figure 6: Fermi L AT abilities to constrain the Milky Way DM density distribution 
after 5 years (keeping correlations with DM particle properties). Same as Fig. Q but for 
Pq = 0.4 GeV/cni^ and = 20 kpc fcvir = 15.2 and M^ir = 1.6 x lO^^ y^^j^ 

which as mentioned above, is the common value used in worlss related to indirect DM detection. 
In this case, Mvir = 1-6 x 10^^ and Cvir = 15.2. Let us note that this value for the scale 
radius implies a concentration parameter much larger than what is obtained, for Mvir ~ 10^^ Mq, 
for the mean in numerical simulations, Cvir ~ 9—11, and it is even away from the lex scatter. 
Thus, although certainly plausible, its choice is not really justified by the value of the mean of the 
concentration parameter obtained in numerical simulations. It is interesting to notice, though, that 
these parameters are very close to those obtained in a very recent estimate |132| . In the left panel 
of Fig. El we see that gamma-ray data alone would not be able to offer any information on the 
DM density distribution in this case either. The addition of external constraints would allow to 
substantially restrict the allowed parameter space, as can be seen from the right panel. However, 
only limited information could be extracted, even in the case of our most optimistic priors. At 90% 
CL (2 dof), the allowed parameter space would be compatible with the mean obtained in numerical 
simulations. 

On the other hand, in Fig. [7]we use a recent estimate of the dark halo mass within the innermost 
60 kpc of the Milky Way |128) . which predicts a lower mass than other works, and we assume 
Mdh(60) = 4.0 X 10^^ Mq. By considering a local DM density given by our default value, /3q = 
0.4 GeV/cm^, this implies a scale radius much smaller, = 10 kpc, than what is usually assumed 
and a concentration parameter correspondingly much larg er, Cvir — 24.4, than what is predicted 
from N-body simulations for halos with virial masses about 10^^ M© (cf. Eq. ([5])). From the left 
panel, we see that a large part of the la scatter band obtained by numerical simulations would 

'^In the case of pQ — 0.3 GeV/cm^ and Ts — 20 kpc, one obtains Cvir = 13.6, which is still much larger than the 
predictions of the mean value in numerical simulations. 
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Figure 7: Fermi L AT abilities to constrain the Milky Way DM density distribution 
after 5 years (keeping correlations with DM particle properties). Same as Fig. Q hut for 
Pq = 0.4 GeV/cni^ and M^h = 4.0 x lO^^ Mq (vg = 10 kpc, c^ir = 24.4 and M^ir = 8.4 x 10^^ Mq). 

be disfavored by such a measurement with only gamma-rays. This trend is further accentuated if 
external information on Pq and (av) is available. In that case, even with our less optimistic priors, 
predictions of numerical simulations would be disfavored at more than 90% CL, as can be seen from 
the right panel. Hence, in this situation a Fermi-hAT detection of DM would point out that a 
revision of our current understanding of halo clustering and evolution is mandatory, perhaps due to 
the influence of baryons, which are usually not included in the simulations. 

All in all, we see that numerical simulations tend to predict a mean concentration parameter 
which typically implies low statistics in a gamma-ray experiment and thus, extracting information 
about the DM density distribution would become a very challenging task in this case. In contrast, 
typical values for the DM local density imply large values of the concentration parameter, which 
decrease with the virial mass. Hence, in order to avoid strong conflict with numerical simulations, 
a virial mass in the high-mass side of the allowed interval would be needed. However, a recent 
determination of the DM halo mass within the innermost 60 kpc |128| points to a low virial mass, 
which would imply a large concentration parameter, in clear conflict with numerical simulations 
(for commonly inferred values of the DM local density), but with much more promising prospects 
to constrain the halo density distribution with a gamma-ray experiment like Fermi-hAT . 

6 Conclusions 

During the last two decades, the progress in high resolution N-body simulations has led to the dis- 
covery of a universal density profile for CDM halos |17pi8|fTTT1 l33 pil5V[118| that can be parametrized 
in analytical forms containing very few parameters. However, these prescriptions only represent the 
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mean results for all simulated lialos, but cannot describe the large scatter that is found |35fl40| . 
Moreover, these simulations do not usually include the effect of baryons, which could significantly 
change the picture |167H175] . Hence, for the particular case of our Milky Way, these mean values 
do not necessarily describe its CDM halo. 

This is a very important issue, as the predictions for the potential signals for the astrophysical 
detection of DM critically depend on the structure of the Milky Way halo. Indeed, several approaches 
have been proposed to determine DM particle properties by using future indirect WIMP DM- 
induced signals of gamma-rays or neutrinos, direct detection searches, collider information or their 
combination |57H108) . All these measurements would be complementary and would constitute an 
important step toward the identification of the particle nature of DM. However, without a better 
knowledge of the Milky Way DM density profile, large uncertainties would be present. This, in 
the case of a positive signal, would make the task of constraining DM properties a much more 
challenging one. 

In the case of gamma-ray studies, although different profiles have been considered in the litera- 
ture, a single DM density profile is commonly assumed to determine the capabilities of gamma-ray 
experiments to extract some DM properties |57H62| . and thus in general, the obtained results are 
too optimistic. Here, we have studied the effect these astrophysical uncertainties would have on the 
determination of some DM particle properties, as annihilation cross section, mass and branching 
ratio into dominant annihilation channels. Conversely, gamma-ray searches could also be used to 
learn about the structure of the Milky Way DM halo, as a complementary tool to other type of 
observational data that study the gravitational effect caused by the presence of DM. 

In this work, we have studied the capabilities of the Fermi-hAT instrument on board of the 
Fermi mission |53) to tackle these issues and consider the potential gamma-ray signal in a squared 
region with a side of 20° around the GC ([/[, |6| < 10°) after 5 years of data taking. In order to 
model the relevant gamma-ray backgrounds, we use the numerical code GALPROP [54] and the 
latest Fermi-L AT observations |55p56). 

In two introductory sections (Sections [2] and [3]) we review the main components of the gamma- 
ray emission from DM annihilation in the GC, and critically discuss the current knowledge of the 
parameters involved in the description of the (smooth) Milky DM density profile. Then, in Section U] 
we describe our modeling of the Fermi-hAT experiment and the ingredients of the analysis we 
perform. Finally, in Section [5l we present the results of our paper. First, we consider an Einasto 
profile and study the Fermi-hAT capabilities to simultaneously constrain DM particle properties and 
the parameters that describe the DM density profile. We do so for two DM masses, = 80 GeV 
(Fig. [T]) and ruy. = 25 GeV (Fig. [5]). These results could be compared to those of Ref. |61| . where 
a particular DM density profile was assumed to perform the fits, but with a single angular bin. In 
addition, we also consider the improvement in these results when external information on (av) and 
Pq is included (Figs. [3]and[l|). In the last part of this work, we focus on the determination of the 
parameters that describe the DM density profile. In order to compare our results to the relations 
obtained in numerical simulations, we consider a NFW profile, study different sets of values for the 
parameters and plot the results in the (Mvir,Cvir) plane (Figs. [SHZ])- The parameters and relevant 
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information for each of the figures are summarized in Table [T] 

In summary, and baring in mind the difficulties of all experiments that aim to detect DM 
to distinguish that signal from any other possible backgrounds, our study shows the Fermi-h AT 
capabilities to simultaneously constrain DM particle properties and the Milky Way DM density 
distribution. Along these lines, we have tried to point out some important issues that should be 
taken into account in indirect searches when a potential DM signal is detected. 
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